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Abstract 

The evolution with carrier concentration of the electrical properties of amorphous indium oxide (InO) thin films has been 
studied using electronic double layer transistor configurations. Carrier variations of up to 7 x 10^** carriers/cm^ were achieved 
using an ionic liquid as a gate dielectric. The superconductor-insulator transition was traversed and the magnitude and 
position of the large magnetoresistance peak found in the insulating regime were modified. The systematic variation of the 
magnetoresistance peak with charge concentration was found to be qualitatively consistent with a simulation based on a model 
involving granularity. 



The competition between disorder and superconduc- 
tivity has long been the subject of theoretical and ex- 
perimental study[l-6]. By varying the nominal disor- 
der a superconductor-insulator (SI) transition is found 
in materials such as amorphous Bi ('a-Bi)[7'|, TiNxQ, 
and InOx[9, 10]. Also Parendo et al. electrostatically 
tuned the SI transition of a-Bi using SrTiOa as both a 
substrate and a gate dielectric [ll'|. Here we focus on the 
electrostatic tuning of amorphous InOx which we will re- 
fer to as InO. This material is of interest because of the 
novel feature of a giant magnetoresistance (MR) peak 
at low temperatures in the insulating re gim e, as first re- 
ported by Paalanen, Hebard, and Ruel[l^, and subse- 
quently studied by many workersflo',n^. Direct evidence 
of localized Cooper pairs, both above the transition tem- 
perature, and in the insulating regime has been reported 
in scanning tunneling microscope (STM)[l3| and planar 
tunneling junction experiments [15|, where the existence 
of a superconducting gap was inferred from the form of 
the density of states. Recent theoretical works (l6l-[l8j 
suggest that the non-monotonic behavior of the MR of 
disordered systems in the insulating state may be associ- 
ated with a competition between the single-grain charg- 
ing energy Ec and the electron paring energy A, modeling 
the film as a granular superconductor. Despite many ex- 
perimental and theoretical works, little is known about 
the role of carrier density in the SI transition of InO and 
the effect of carrier density variations on the MR peak. 
Here we present data relating to these issues and sug- 
gest a qualitative explanation of the results based on a 
numerical simulation using one of the models(l8|. 

Amorphous InO films were grown on 0.6 x 0.6 cm^ 
SiO substrates and were patterned during growth using 
a physical (shadow) mask. Sample thicknesses ranged 
from 10 nm to 6 nm. The 6 nm samples were all insulat- 
ing as grown and are the central focus of this work. Gold 
electrodes, 10 nm thick, were deposited on the InO films. 
An ionic liquid (IL), DEMI3-TFSI, was used as a gate di- 
electric with a Pt coil top gate to form an electric double 



layer transistor (EDLT). A schematic of the top-gated 
EDLT can be found in our previous work[19]. (Earlier, 
Misra, McCarthy, and Hebard gated InO films with ILs 
but did not explore their superconductivity [2(]| .) 

The sheet resistances Rs of the films were determined 
employing a four-probe electrode configuration. A ■^He 
refrigerator with a superconducting magnet enabled us 
to vary the temperature T between 300K and 0.40K and 
the magnetic field H between OT and 9T. Gate voltages 
Vg, ranging from +3V to -3V were used to induce or 
deplete charge carriers, which are electrons. The various 
values of Vg were applied at a temperature of 240K and 
were held constant throughout the subsequent cooling 
and measurement. 

The T dependencies oiRs, of two InO samples at differ- 
ent values of are plotted in Fig. 1. Both samples were 
grown under the same conditions i.e., at the same oxy- 
gen pressure and at the same time. However their initial 
resistances were adjusted so as to be different, by anneal- 
ing at 65°C under vacuum for different periods of time. 
Both samples exhibited large changes in Rg upon gating. 
Unfortunately, Hall effect measurements of such highly 
disordered systems were not possible, except in the case 
of a few other samples (not shown) that had relatively 
low normal state resistances. The process of gating was 
found to be reversible, with minor hysteresis. This sug- 
gests that apart from the remote possibility of a reversible 
chemical reaction, the gating process is electrostatic. To 
determine and/or confirm the carrier modulation of the 
system, we used an electrochemical technique known as 
the Chronocoulometryjilj. Using this technique, we ob- 
served changes in the sheet (2D) carrier density or charge 
transfer (Aris/ieet) of up to 7 x IQ^"^ carrier s / cm? . How- 
ever we did not characterize the charge transfer at each 
gate voltage for the films reported here. 

The accessible charge transfer was large enough for us 
to observe the SI transition of InO as shown in Figs. 
1(a) and 1(b). The right hand panels of Fig. 1 show Rg 
vs. T behavior of the insulating state; both Arrhenius 
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Figure 1: (Color online) The temperature dependence of Rs 
at various Vg of samples A (a) and B (b). The onset of su- 
perconductivity (the initial downward dip in the Rs vs. T 
curve) occurs at 0.85K, with Vg = +1.001/ for sample A, and 
at 1.47K, with Vg — —2.00V for sample B respectively. On 
the right hand panels, both Rs vs. 1/T and l/T^^^ are plot- 
ted for sample A (al-a2) and for sample B (bl-b2). Note the 



cross-over from Arrhenius to ES VRH at Vg 



-0.50y shown 



in (al-a2), whereas (bl-b2) shows ES VRH over all ranges of 

Vg. 
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Figure 2: (Color online) Rs vs. H at various values of T and 
Vg for sample A, labeled (al-a3) and sample B, labeled (bl- 
b3). The density of electrons increases from the left to the 
right. Each curve represents a separate isotherm of Rs vs. 
H, where T ranges from 0.5~1.3K for sample A (more disor- 
dered, top) and from 0.5~0.9K for sample B (less disordered, 
bottom). 



R = RoexpiTo/T) and Efros-Shklovskii variable range 
hopping (ES VRH) R = Roexp{To/T)^/'^ were observed. 
Comparing Figs. l(al) and l(bl), the more insulat- 
ing sample shows a wider range of Arrhenius activated 
transport, suggesting that it has a hard gap in its single- 
particle density of states. A cross-over from Arrhenius 
to ES hopping is observed as the film becomes less resis- 
tive with increasing carrier concentration. The observed 
ES hopping suggests the presence of a soft gap (Coulomb 
gap) in the density of states near the Fermi level, which 
is due to long range interactions between electrons in the 
system (i^. This has also been reported in other granu- 
lar systems[^ [2^. Both samples A and B exhibited SI 
transitions tuned by carrier modulation, which is clearer 



in the case of sample B [Fig. 1(b)]. 

The MR measurements of these two samples at differ- 
ent gate voltages are shown in Fig. 2. As mentioned, such 
MR peaks in disordered systems have been reported in 
a variety of experiments[3. [lol. [l^. [l3| . The most signifi- 
cant feature of the data for sample A [Figs. 2(al), 2(a2), 
and 2(a3)] is the transition from negative MR to positive 
MR followed by downward slope of Rg (MR peak) upon 
further increase of , all taken at fixed T over the range 
from 0.5K to IK. 

We have attempted to quantify some of our observa- 
tions of the MR peak. A magnetic length ( associated 
with the field at the peak Hpeak, can be estimated as 
( = ($o/ffpeafc)^/^ where $o = 2.07 x 10"^^ T-m^ is the 
flux quantum. Values of the peak magnetic field, Hpeak 
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at different gate voltages, all taken from the isotherm 
curves of Rs vs. H aXT = Q.'bK , were found to be mono- 
tonically increasing functions of the gate voltage. This 
implies that C decreases with increasing carrier density. 
A measure of the size of the MR peak can be taken to be 
Rnorm = Rs{Hpeak) / Ra{H = 2Hpeak)- Figure 3 shows a 
plot of Rnorm'vs. ^ for cach of the samples, where differ- 
ent points correspond to different values of gate voltage. 
The inset to Fig. 3 shows a plot of Hpeak vs. the recipro- 
cal of Rs at lOOK. The latter is proportional to the carrier 
concentration assuming a Drude resistance at high tem- 
peratures and assuming that the gate voltage changes the 
carrier concentration and not the disorder [2J|. Note that 
in the inset, curves of Hpeak vs. 1/Rs,iook of the two 
samples collapse onto each other, suggesting a possible 
universal behavior. 

There is an apparent connection of this behavior to 
sample morphology. Atomic force microscope (AFM) 
images (Fig. 3, bottom, left) of these films exhibited 
roughness which could be interpreted as granularity, even 
though the films are completely connected. A Power 
Spectral Density (PSD) analysis (Fig. 3, bottom, right) 
of the AFM images revealed a structural length scale 
implying a weak level of correlation in the morphology. 
These lengths are indicated by the vertical dotted lines 
with estimated error bars for the two curves of Fig. 3. 
Note that the lengths corresponding to the PSD peaks 
are close to the values of C corresponding to the maxima 
of Rnorm'vs. ( especially in the case of sample A. Thus 
when the magnetic length scale corresponds to the struc- 
tural length, there is a maximum in the strength of the 
MR peak. 

A granular morphology may not be essential for a sys- 
tem to exhibit granular like behavior. The oxygen con- 
centration may be the key to determining an effective 
granularity in the amorphous InO system. A recent ex- 
perimental workf^s'l showed that mesoscale spatial vari- 
ations in the oxygen concentration are present in nomi- 
nally homogeneous, amorphous InO films. Such chemical 
inhomogeneity can induce local carrier density fluctua- 
tions resulting in mesoscale inhomogeneity of the pair- 
ing energy A. However, we don't know if the mesoscale 
thickness variations in the films of the present work are 
correlated with variations in the local indium to oxygen 
ratio. 

In a disorder-driven SI transition, the spatial inhomo- 
geneity of the pairing energy A can be very important. 
In the insulating regime, in some models, the system 
may break up into superconducting islands, i.e., become 
effectively granular [^. In this picture global supercon- 
ductivity occurs when there is phase coherence between 
the localized order parameters of different grains, while 
the energy to break a pair is much larger and remains 
nonzero even when the coherence (and thus the global 
superconductivity) is destroyedfT^. In Ref. \d\ such a 
spatial fluctuation of the order parameter amplitude and 
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Figure 3: (Color online) (Top) The plot of R„ormVS. ( for 
samples A(red, square) and B(blue, triangle). The inset shows 
the values of peak field Hpeak where the maximum of the MR 
at T = 0.5K is found as a function of l/-Rs,ioox in zero mag- 
netic field, which is proportional to the carrier concentration. 
The quantity defined as Rnorm is a measure of the magnitude 
of the MR peak. The maximum Rnorm occurs very near the 
length extracted from the PSD of the AFM surface profile 
(dotted vertical lines) for sample A and very close to that 
length for sample B. (Bottom) The AFM image (left) and 
the PSD analysis (right) showing a peak, corresponding to 
the periodicity of surface morphology, for sample B. Similar 
results were observed for sample A. 



corresponding phase correlation have been demonstrated 
using a negative U-Hubbard model. Near the SI tran- 
sition, it has been suggested [26* 



27| that the MR peaks 



in disordered systems arise because magnetic flelds affect 
the concentration and size of superconducting islands, so 
that as these islands shrink with increasing held there is a 
transition from Cooper pair-dominated to single electron- 
dominated transport. 

On the other hand, reduction of the superconduct- 
ing pairing energy A within grains can itself lead to a 
tradeoff between conduction by Cooper pairs and con- 
duction by unpaired electrons, and thus (potentially) to 
a MR peak, even when the concentration and the size of 
the superconducting islands are fixed. Recent theoretical 



works [171, Il8j study a model with fixed size and concen- 
tration of superconducting grains, presumably set by the 
film morphology, and show how a MR peak deep in the 
insulating state can arise as a result of the reduction of 
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Figure 4: (Color online) Simulation of the resistance of a 2D 
array of identical superconducting grains deep in the insulat- 
ing state as a function of H. Different curves are labeled with 
their corresponding values of Aq/Ec, and are normalized to 
their maximum Rmax- As Ao is increased, which presumably 
corresponds to larger carrier density, a MR peak develops that 
shifts to larger magnetic field, in qualitative agreement with 
what is seen in Fig. 2. Here all curves correspond to a temper- 
ature such that fcsT = O.lEc and have localization lengths 
and ^2 for single-electron and pair conductivity, respectively, 
satisfying ~ 8. 

the superconducting gap A with increased H. This pre- 
dicts that near the MR peak and at low temperature the 
conduction should be described by ES VRH, as shown in 
Figs. I(a2) and l(b2). Both approaches lead to an insu- 
lator in which Cooper pairs with nonzero A are formed in 
the insulating regime of the system and are responsible 
for the MR peak. 

The shift of the MR peak to higher magnetic fields with 
increasing carrier concentration, as shown in Figs. 2 and 
3, can be explained qual itatively within the context of 
the theory of Ref . |l6l4lg | . Increasing the carrier density 
presumably increases the density of states at the Fermi 
level within the superconducting grains, thereby driving 
up the zero-field superconducting gap, Aq. A larger Aq 
implies that a larger H is required in order to reduce A 
to the value of the grain charging energy Ec, so that the 
MR peak shifts to higher H. In this way the transition 
from negative MR [as in Fig. 2(al)] to a peak at an 
intermediate H [Fig. 2(a2)] to a peak at a large H [Fig. 
2(a3)[ can be understood. 

As an example. Fig. 4 shows values of the resistance of 
a simulated 2D array of regularly-spaced, monodispersed 
superconducting grains as a function of 7J , calculated us- 
ing the method described in Ref. At small Aq/ Ec, 
the conductivity is primarily due to hopping of unpaired 
electrons, and there is a monotonic negative MR [as seen, 
for example, in Fig. 2(al)[. At larger Aq/Ec, which os- 
tensibly corresponds to larger carrier density, the MR 
develops a peak associated with a trade-off between con- 
ductivity by single electrons and conductivity by Cooper 
pairs. This peak moves to larger H as Aq/Ec is increased 



[as in Fig. 2]. For the simulation of Fig. 4 we have as- 
sumed a conventional BCS-like dependence of A on the 
field H: A ^ Aq ^/T~{H/Hc)^ . In this way the data 
shown in Fig. 2 is consistent with the concept of tuning 
the local superconducting gap by modulating the carrier 
density. Unfortunately, the simulation method used to 
generate Fig. 4 cannot be used for a quantitative de- 
termination of the relationship A(n), since this requires 
a knowledge of the iJ-dependence of the gap as well as 
the relative localization lengths ^1 and ^2 for unpaired 
and paired electron hopping. We also caution that the 
simulation technique is applicable only for the heavily- 
insulating limit, and in this sense our comparison be- 
tween Figs. 2 and 4 is only qualitative. It should also be 
noted that within this simple model a strong MR peak 
develops only at relatively large £,2/^1- A final caveat is 
the possibility that other models may give similar results. 

In conclusion, we have demonstrated the evolution of 
the MR peak found in InO as a function of carrier density. 
The transition from a negative MR to a strong MR peak 
followed by the suppression of the MR peak was observed 
near the SI transition tuned by carrier modulation. Our 
findings support the presence of localized Cooper pairs 
in the insulating regime and are qualitatively consistent 
with models based on granularity of the system. The 
reason the amplitude of the MR peak is a maximum when 
the magnetic length is close to a length scale associated 
with disorder in the films remains an open question. 
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